function x = UnitUxbSol(U,b,n)
%UnitLxbSol 解单位上三角阵方程
%输出结果x储存在b中
for i = n:-1:2
    b(1:i-1) = b(1:i-1) - b(i)*U(1:i-1,i);
end
x = b;
end

